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ABSTRACT 

We investigate the impact of instrumental systematic errors in interferometric measurements of 
the cosmic microwave background (CMB) temperature and polarization power spectra. We simulate 
interferometric CMB observations to generate mock visibilities and estimate power spectra using the 
statistically optimal maximum likelihood technique. We define a quadratic error measure to determine 
allowable levels of systematic error that do not induce power spectrum errors beyond a given tolerance. 
As an example, in this study we focus on differential pointing errors. The effects of other systematics 
can be simulated by this pipeline in a straightforward manner. We find that, in order to accurately 
recover the underlying B-modes for r = 0.01 at 28 < £ < 384, Gaussian-distributed pointing errors 
must be controlled to 0.7° rms for an interferometer with an antenna configuration similar to QUBIC, 
in agreement with analytical estimates. Only the statistical uncertainty for 28 < I < 88 would 
be changed at ~ 10% level. We also show that the impact of pointing errors on the TB and EB 
measurements is negligibly small. 

Subject headings: cosmology:observations, cosmic microwave background, instrumenta- 
tiominterferometers, methods: data analysis, statistical 



1. INTRODUCTION 

Cosmic microwave background (CMB) polarization 
measurements can significantly improve the estima- 
tion of cosmological parameters, breaking the degen- 
eracies between parameters measured using CMB tem- 
perature anisotropy data alone. In the standard the- 
ory of the CMB, the polarization field can be de- 
composed uniquely into an electric-type E-mode and 
a magnetic-type i?-mode (|Zaldarriaga fc Seliakl Il997t 
iKamionkowski et al.l Il997ah . The E-mode polarization 
can provi de useful informatio n about reionization of the 
universe (iHu fc Holder! [20031 ) . The primordial _B-modes 
can probe horizon-scale primordial gravitational waves 
and play a major role in unders tanding the inflation- 
ary epoch ()Hu fc Dodelsonl l2002t ) , while the secondary 
lensing-induced E-mode signals (|Zaldarriaga fc Seliakl 
promise to provide a wealth of information about 
the distribution of matter and the evolution of large scale 
structure. Measuring the CMB polarization has become 
one of the major goals of CMB experiments. However, 
the polarized CMB signal is so small that its measure- 
ment requires not only very high instrumental sensitivity, 
but also exquisite control of systematics. 

In many traditional imaging experiments, the deter- 
mination of the Stokes parameters Q and U is based on 
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subtracting intensities measured by two different detec- 
tors. S uch an experime nt is very sensitive to systematic 
errors (|Hu et al.l |2003|) . For instance, beam imperfec- 
tions or beam mismatch will cause leakage from total 
intensity I into polarization s ignals Q and U. A re- 
cent study ijMiller et al.l 120081) shows that in order to 
achieve a reliable E-mode detection (r = 0.01), allow- 
able levels of beam systematics should not exceed 1% 
in ellipticity, the sub-percent level in differential beam 
width and the few- to sub-arc sec level in differential 
pointing. Also, with a finite patch of sky observed by 
single-dish intruments, it is impossible to perfectly sep- 
arate the very weak E-modes from the m uch stronger 
E-mo des dlewis et al.l [20021: iBunnl |2002al lbl: iBunn et all 
120031 : IBunnl 1201 ll ). Therefore polarization detection 
presents a great challenge in imaging experiments. 

Alternatively, interferometers are a more natural 
choice for measuring the anisotropics of the CMB tem- 
perature and polarization. The correlation of the electric 
fields from two antennas, called a visibility, measures the 
Fourier transform of the sky intensity fluctuations modu- 
lated by the response of the antennas. In most cases the 
region of sky covered by the antennas is small enough 
that one can use the "flat-sky" approximation. The ex- 
pansion of the intensity field into spherical harmonics 
thus can be approximated by Fourier modes - the visi- 
bility directly relates to the CMB power spectrum. The 
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main reason for building interferometers instead of tra- 
ditional imaging experiments is that systematic effects 
are well-co ntrolled in s ome cases, especially for B-mode 
detection (jBunnl 120071 ). Since interferometers measure 
the Stokes parameters directly - without differencing the 
signals between separate detectors for measuring polar- 
izations - mismatched beam shapes and pointing errors 
do not cause leakage from / into Q and U. In addition, 
in contrast to imaging experiments, interferometers can 
separate the E- and B-modes more cle anly because they 
sample the sky in the Fourier domain ([Park et alj|2003t 
iPark fc Ngll200l . 

Interferometers have been used to measure the CMB 
anisotropics since the 1980s. The first attempt at 
measuring t he CMB with an in terferometer was car- 
ried out by iMartin et ail (|1980[) . Shortly after that, 
the 27 anten na VLA was used in searches for CMB 
fluctuations (iFomalont et all Il984t iKnoke et all 11984 
iPartridge et al.lll988l ) and the first dedicated interferom- 
eter with two-eleme nt correlation receiver to C MB re- 
search was made bv lTimbie fc Wilkinson! (|1988l ). So far 
a number of interferometers have been constructed to 
observe the CMB power spectrum. The CAT telescope 
was the first interfe rometer to actually detect struc- 
tures in the CMB (lO'Sullivan et all 119951 IScott et all 
119961: iBaker et all 1 19991) and DASI (jKovac et all 12002]) 
first detected the faint po larized signa ls in the CMB. 
CBI dPearson et al l 120031) and VSA (jDickinson et all 
l2004tlGrainge et al.ll2003| ) made high-sensitivity observa- 
tions of the CMB temperature and polarization angular 
power spectra down to sub-degree scales. 

There are many papers in the literature on the study 
of how instrumental systematic errors affect CMB an- 
gular power spectrum measurements. For imaging mea- 
surements, a pioneer i ng stu dy of such effects was per- 
formed by iHu et all (|2003f ). In addition, contamina- 
tion of the CMB power spe ctrum by systematic effects 
has been precisely assessed (ISu et al 1201 lb iMiller et aT 
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20101) . An analytic approach for characterizing a vari- 



ety of syst ematic errors for interferometers has been per- 
formed by iBunnl (J2007) . This approach diagnoses sys- 
tematic errors in a qualitative way. Actual experiments 
will naturally require a realistic simulation to quantita- 
tively study such effects as carefully as possible. In this 
paper, we present for the first time a simulation pipeline 
to accurately assess the impact of interferometric sys- 
tcmatics on CMB power spectrum measurements. Such 
systematic errors are evaluated from a full maximum like- 
lihood (ML) analysis of realistic simulated data. The 
method presented in this paper is able to characterize 
a wide variety of systematic errors, such as beam shape 
errors, gain errors and cross-polarization, etc. As an ex- 
ample, we examine a specific configuration of an inter- 
ferometer to quantify systematic pointing errors on the 
recovery of the CMB power spectrum. 

The remainder of this paper is organized as follows. 
In Sec. [H we briefly summarize the simulation of CMB 
interferometric observations. In Sec. we describe the 
ML method for extracting the CMB temperature and 
polarization power spectra from interferometric data. In 
Sec. 01 we focus on pointing errors. Using a comparison 
of the recovered CMB power spectra with and without 
pointing errors, we assess the degree of contamination 



of such systcmatics on CMB power spectrum recovery. 
Finally, a discussion and summary are given in Sec. [5j 

2. SIMULATIONS 
2.1. Visibilities and Covariance Matrix 



Following the p revio u s papers (|Bunnl 1 2007 



: Maismger \2UUJt Iriopson & Maeuend UUMbt 
Mvers et all 120031 . 120061: IPark et all 120031 : iWhite eTall 
19991 ). here wc briefly review CMB interferometric 
observations. Suppose that we have a set of antennas 
at position r n , n = 1, . . . , and each antenna measures 
two polarization states (linear polarizations or circular 
polarizations). The signal received, e ou t, by the n-th 
antenna in response to incoming radiation fields, £.;„, at 
frequency, v, from direction k is then 



d 2 kG(k)-e m (k)e 



i(k-r n —1-Kvt) 



(i) 



Here G is the 2x2 matrix-valued antenna pattern 
which encodes the antennas' response to the sky. £;„ 
and e ou t are two-dimensional complex vectors represent- 
ing an incoming polarized electric field and an electric 
field that is output by the antenna, respectively. The 
two-component vectors e\ n and t out can be expressed 
in either a linear polarization basis (X-Y) or a circular- 
polarization basis (R-L). These two bases are connected 
by a unitary transformation. 

The time-averaged value of the correlation between po- 
larization component m from antenna j, and polarization 
component n from antenna k is referred to as a visibility: 

Le - Vmn = ( e lut, m e out,n) ■ In the flat-sky approximation, 
for both the linear and circular-polarization bases, the 
visibilities with a 2 x 2 matrix form are related to the 
Stokes parameter matrix as follows: 

V 3k = J d 2 xG J (x) • S(x) • G+ fe (x)e- 27Hu '"- x , (2) 

where the baseline vector Ujk measures the separation 
between the two antennas in units of wavelength (rj — 
rfc)/A. The Stokes matrix, S, is related to the Stokes 
parameters I,Q,U,V for linear and circular-polarization 
bases: 

' I + Q U + iV 
U-iV I-Q 



(3) 



Scire — 



I + V Q + iU 
Q-iU I-V 



(4) 



Us ually two m ain descriptors for systcmatics are 
used (|Bunnll2007l ) : the instrumental Jones matrix, J, and 
the antenna pattern, G(x). J describes systematics that 
are purely introduced within the instrument, such as gain 
errors and cross-talk between the two outputs of a given 
antenna. Instrumental errors are easier to model since 
they do not depend on the position on the sky, while 
G(x) characterizes systematics occurring from the ob- 
servation of the sky before instrumental errors are taken 
into account. We thus use G(x) to model beam shape 
errors, pointing errors and cross polarizations, etc. This 
separation is a bit arbitrary in that instrumental errors 
could be absorbed into antenna patterns, but they are 
a convenient conceptual distinction between systemat- 
ics happening "before" and "after" the antenna averages 
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over the beam. The total effect on the visibilities can be 
found from the relations: 

V lk = J d 2 x J 1 (G J (x) ■ S(x) ■ G tfe (x)) jtk e -2*iu Jfc .x . 

( 5 ) 

For the j-th antenna with gain errors g\ , g\ and cross- 
talk couplings e\, e J 2 , the Jones matrix reads 



1 



■A 



i 



si 



(6) 



Similar to the Jones matrix, an azimuthally symmetric 
antenna pattern has the form of 



G J (r,< 



iGjcos(20) ±Gjsin(2<; 



±Gj sin(20) 



G£ - \G{ cos(2(/>) 



(7) 

where (r, </>) are polar coordinates, Go is the ideal beam 
shape and G\ leads to two polarization states mixing. 
The scalar functions Go, Gi depend only on r. For 
an ideal interferometer, each antenna has an identical 
response to both polarization states while there is no 
mixing between them. In this case, G is equal to a 
scalar function multiplied by the identity matrix, i.e., 
G(x) = G(x)l. 

For both linear and circular polarization experiments, 
coupling errors (e) are the major sources of systemat- 
ics affecting the B-mode power spectrum due to leakage 
from / into Q and U, while less worrisome gain errors 
would mi x Vq and Vjj with each other. Furthermore, as 
shown by iBunnl (|2007t ) . the visibility for Stokes V could 
provide a useful diagnostic to monitor the presence of 
these systematics. 

For simplicity we assume that instrument errors in the 
Jones matrix are negligible (taking J = 1) and each an- 
tenna has identical beam patterns for both polarization 
states and has no cross-polar response (i.e., off-diagonal 
entries in G). Then, each visibility measures a simple 
linear combination of the Stokes parameters. We can 
extract Stokes visibilities from Eq. [2J yielding 



V 3 z k = J d 2 zZ(x)G J (x)G* fc (xK 



(8) 



for Z = {I,Q,U,V}. Since the visibility function is the 
Fourier transform of the Stokes fields on the sky weighted 
by the antenna response, by using the well-known con- 
volution theorem, we can write the visibility function in 
Fourier space (wv-domain) : 



(9) 



where Z(u) and G 3 (u) are the Fourier transforms of 
Z(x) and G- 7 (x)G* fe ( x ), respectively. Note that each an- 
tenna pattern could differ from the others because of sys- 
tematic beam errors. 

In this study, we only focus on differential pointing er- 
rors. The pipeline presented in this paper can be applied 
in a straightforward manner to other systematics. Point- 
ing errors occur when not all the anten nas point in the 
same direction. Following IBunnl (|2007t) , we model the 
pointing error as a Gaussian-distributed error with dis- 
persion S. Assuming the j-th antenna has the pointing 



offset (5xj relative to a desired direction, then according 
to Eq. [5J each visibility can be expressed in the form 



V, 



J d 2 xZ(x)G 1 (x + SiL J )G* k (x + 5x k ) 



V z (u) = J d 2 xZ(x.)A(x.)e- 2 ™ 
Vz(u) = I d 2 w Z(u - w)l(w) 



(10) 

where Z = {I,Q,U} and the beam response G(x) can 
be approximated by a circular Gaussian of a frequency- 
dependent dispersion tr, i.e., G(x) = exp(— |x| 2 /2<r 2 (^)). 
The corresponding Fourier transform is thus given by 
G(u) = 27R7 2 (i/)exp(-27r 2 |u| 2 CT 2 H). 

In order to recovery the power spectrum based on sim- 
ulated visibility data, one needs to construct the covari- 
ance matrix, which is the fundamental tool for analysis 
of Gaussian random CMB fields. In principle, all kinds of 
systematic errors in visibilities can be simulated through 
Eq.[5l whereas the theoretically predicted covariance ma- 
trix does not include any systematic uncertainties. Tak- 
ing the error-free beam pattern, G, with identical re- 
sponse between antennas, Eqs. [5] and HO can be simplified 
as follows: 



(11) 
(12) 



Here the intensity beam pattern, A(x), is defined by 
|G(x)| 2 and A is its Fourier transform. 

In the flat-sky approximation, the Stokes parameters 
Q and U can be decomposed into the E- a nd £?-modcs in 
Fourier space (jZaldarriaea fc Seliak|[T997t) . Using Eq.[H 
the Stokes visibilities Vj , Vq and Vjj then can be ex- 
pressed in terms of T, E and B modes as follows: 

V Q (u)= J d 2 w[E(w) cos(2</) w ) - B(w) sin(2<£ w )]I(u - w) 
Vu(u)= [ d 2 «;[i?(w)sin(20 w ) + i?(w)]cos(20 w )i(u- w) 

(13) 

where T, E and B stand for the CMB temperature field, 
T, and polarization fields, E and B, in Fourier space and 
W is the angle made by the vector w with respect to the 
x-axis. In this study we assume the Stokes visibility Vy 
to be zero. 

Given a set of visibility measurements, one can use 
maximum likelihood analysis to evaluate the CMB power 
spectra of the temperature and polarization. By defin- 
ing a vector of data V = (Vj 1 , V%, Vj); • • ■ ; V?,VZj, V$) 
at each baseline vector with i = 1, ...,n, the corre- 
sponding covariance matrices of the CMB visibilities are 



Vi(u)= / d 2 wf(w)A(u-w), 



d 2 w / d 2 w' (Z(xv)Z~'*(w') 



x A(ui - w)A*(uj - w') 
= J d 2 W S^(|w|)i(u, - w)i*(uj - w) , (14) 

where i,j denote visibility data indices and the depen- 
dence of the correlation functions S zz > on the CMB 
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power spectra are listed in Tabic [T] In the flat sky ap- 
proximation, the 2D powe r spectrum 47r 2 |u| 2 iS'(|u|) ~ 
£{t + l)C^ =27ru for £ > 10 (|White et al.lfl999l) . 

2.2. Simulated Observations 

The CMB Stokes fields are believed to be is otropic 
and Gaussian in the standard infl at ionary models (iGuth 



1981 tlKamionkowski et al.lll997al lb1: rZaldarriaga fc Seliak 



1997|). On a small patch of the sky, the correspond- 
ing Fourier components of these fields are complex ran- 
dom variables and the value of the real and imaginary 
parts of each point u in Fourier space are drawn inde- 
pendently from a normal distribution with zero mean 
and variance cx C^|^ =2w | u |. With cosmological parame- 



ters deriyed_lrom_WMAP_^7-year results ([Larson et al.l 
[20111: iKomatsu et all |20H . ' we use the public code 
CAMB ([Lewis et al.ll2000l ) to compute the CMB power 
spectra Cj T ,Cf E ,Cj E : Cf B . The input B-modc con- 
tains a primordial component with a tensor-to-scalar per- 
turbation ratio, r, and a secondary component induced 
by lensing. In our simulation we fix r = 0.01, the goal 
for many current observations. 

Based on these power spectra, we generate Fourier 
modes and then perform the inverse Fourier transform 
to obtain real-space Stokes fields I(x),Q(x) and U(x.). 
From Eq. [5[ for a given Jones matrix and beam response 
G, the Stokes visibilities are then obtained by performing 
the Fourier transform again. 

We assume the instrumental noise at each point 
of the w-plane is a complex, Gaussian-distributed 
numb er which is indep e ndent between different base- 
lines (jWhite et al.l 119991: iMorales fc Wvithd 120101 ) . For 
an instrument which measures both polarizations with an 
identical uncertainty in Stokes parameters, we can sep- 
arately generate the Gaussian noise with identical rms 
levels for each /, Q and U visibility. The correlation 
function of the noise for baselines i and j is determined 

by 



r< l 3 



= ( 



XT, 



Sy,S 



\r) A A D J \AvtaUb 



1 



(15) 



where T sys stands for the system noise temperature, A 
for the observing wavelength, Ajj for the physical area 
of a antenna, t\a for the aperture efficiency, A„ for band- 
width, rib for the number of baselines with the same base- 
line vector u and t a for the integration time of the base- 
line. 

In order to illustrate the effect of systematic errors on 
the recovered CMB power spectra and set allowable tol- 
erance levels for those errors, we perform simulations for 
a specific interferometer design. We choose an antenna 
configuration similar t o that of the QUBIC instrument 
(Battistell i et al.ll2011l ) which is under construction for 
observations at 150 GHz. In our simulation, the interfer- 
ometer is a two-dimensional square close-packed array of 
400 horn antennas with Gaussian beams of width 5° in 
the intensity beam pattern A(x), corresponding to ~ 7.1° 
in G(x). The antennas have uniform physical separations 
of 7.89A. With this configuration, the resolution in the 
uw-plane is about a u = 1.82 (A^ ~ 11), and the uv cov- 
erage reaches down to £ > 50 — 2A£ = 28, probing the 
primordial i?-mode bump at I « 50. 

We also assume that all Stokes visibilities. /, Q and 
U, can be measured simultaneously for each antenna pair 



with an associated rms noise level of 0.015/iK per visi- 
bility, roughly corresponding to low-noise detectors each 
with 150^Ks 1 / 2 and a total integration time of three 
years. With this noise level, the simulations show that 
the averaged overall signal-to-noise ratio (SNR) in Stokes 
Q and U maps is about 5. The high SNR ensures an accu- 
rate recovery of the i?-mode power spectrum and allows 
us to see systematic effects clearly. 

We generate realizations of Stokes parameter maps 
having a physical size of 30 degrees on a side and resolu- 
tion of 64 x 64 pixels. This large patch size ensures that 
the intensity beam pattern |G| 2 at the edges decreases to 
~ 1-percent level of its peak value. Although this size of 
patch seems to severely violate the flat-sky approxima- 
tion, the primary beam pattern itself is small enough (the 
field-of-view fi is about 0.047 sr) so that the flat-sky ap- 
proximation is still valid. For simplicity, we assume that 
all the antennas continuously observe the same sky patch 
at a celestial pole and the interferometer is located at the 
north or south pole , the uu-tracks should be perfectly 
circular for a 12-h observation. Fig. 12.21 shows the mock 
systematics-free visibility data from these observations. 

3. MAXIMUM LIKELIHOOD ANALYSIS 

The maximum likelihood estimator of the power spec- 
trum has many de sirable properties ([Bond et al.l 119981 : 
iKendall et al.l fl987T ). The idea is to choose a model for 
the data and construct a likelihood estimator to evaluate 
how well the model matches the data. For a given model, 
comparing to the actual data set will give a likelihood of 
the model parameters. In practice, it is easier to max- 
imize the logarithm of the likelihood function than the 
likelihood function itself. 

Since CMB Stokes visibilities are complex Gaussian 
random variables with zero mean and dispersion Cy + 
Cat, the logarithm of the likelihood function is given by 

In C(C t ) =n\ogTT- log \C V + C N \ - V* (C v + Cat)" 1 V , 

(16) 

where V is the visibility data vector, Cy is the signal co- 
variance matrix predicted by (V"I"V), which can be con- 
structed through Eq. [14J and Cn is the noise covariancc 
matrix, which can be computed by Eq. 1151 

In practice, we parameterize the CMB power spec- 
trum Cg as flat band-powers over some multipole range 
to evaluate the likelihood function dBunn fc Whitdll996l : 
iBond et al.lfl998l : lGorski"et~aT1ll996l: iWhite et al.lll999¥ . 
We divide the power spectrum £(£ + l)Cg into Nb 
piecewise-constant bins. Each bin corresponds to sep- 
arate annuli in the itu-planc, characterizing the aver- 
aged Cg over its bin-width. In our case, we evaluate 
the likelihood function by varying the CMB band-powers 

{Cp^Cp^Cp,^} with b = l,...,N b . 
Here C b = 27r|u 6 | 2 S(|u b |). 

The bin-width can be chosen arbitrarily, but an ap- 
propriate choice of width is fine enough resolution to 
accurately detect the structure of the power spectrum 
and also wide enough to reduce the correlation between 
the band-power estimates so that the statistical errors 
on different band-power bins are approximately uncor- 
related. The natural choice of bin- width can be ap- 
proximated by the characteristic width of the Fourier 
transformed intensity beam pattern ^4(x), which defines 
the typical correlation length in the itu-planc. The 
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TABLE 1 

Dependence of ensemble-averaged Stokes parameter correlations on the CMB angular power spectra 



ZZ' 


(Z (w)£'*(w')) = 5„/(|w|)5(w - w') 






II 


5/7 = Stt{w) 






IQ 


Siq = S T e{w) cos 20 w - S T b sin20 w 






IV 


S IV = S T e(w) sin 20 w + S T b sin 20 w 






QQ 


Sqq = S E e{w) cos 2 20 w + S SB («)) sin 2 20 w 


- S E b(w) sin40 w 




QU 


>5q[/ = (See{w) - Sss(«j))sin2</i w cos2(/) w - 


h Sbb (cos 2 20 w — sin 2 


20w) 


UU 


5(7fj = -See (tu) sin 2 20 w + Sbb (w) cos 2 20 w 


+ 2Seb (w) sin 20 w cos 


2</-w 



TABLE 2 

Table of integrals used in the calculations of the window functions. 



f 

Jo 



exp(47r 2 q • w 



f n d4>vj exp(47r 2 q ■ w 
/ 27r d0 w exp(47r 2 q • w 
J* 27r d0 w exp(4ir 2 q • w 
J 27r c(0 w exp(47r 2 q ■ w 
J*„ 27r (i</> w exp(47r 2 q • w 



= 2vr/ (a) 

cos20 w = 271/2(1) cos(2</3 q ) 
sin20 w = 2nl2(a) sin(2</3 q ) 
cos 2 20 



sin 2 20 v 



7r(a7 (a)((a 2 +24) cos(40 q ) + a 2 ) -8(a 2 +6)ii (a) cos(4^.q)) 

^ 

7r(8(a 2 +6)'i(a) cos(44>q) + / (a) (a 3 -a(a 2 +24) cos(44>q))) 



sin 20 w cos 20 v 



7r(a( t i 2 + 24)7 (a)-8( [ i 2 +6)/i(ii))ain(4</> q ) 
: 



where we introduce |q| = cr(ui + uj), a = 47r 2 |q||w| and I m (a) is the modified Bessel function 
of the first kind and order m for a real argument a (see details in the text). 



minimum bin-width for the QUBIC-likc experiment is 
about Au « VS! -1 =4.1 wavelengths, corresponding to 
A£ sa 26. As a consequence, the total number of band- 
power bins is 6 x Nb ~ 84. However, the computational 
time required to evaluate the likelihood function in such a 
large number of bins is unfeasible. Instead, in this paper, 
we estimate the power spectrum by using the bin-width 
of A£ ~ 60, roughly having 6 band-power bins for each 
power spectrum at the range of 28 < £ < 384. 

Using the abo ve parametrization and fo l lowing the 
previous papers (iHobson fc Maisingerll200l iPark et al.l 
[20031: iWhite et al.lll999ft . the covariance matrices defined 
in Eq. Q3] can be written as 



r*3 
u zz> 



N, 



b=l a,/3 



a/3 



l u f>ll 



where we introduced the so-called window functions 

w zz> af 3 s iven b y 



d(j) w uj Za uJz'i3M u i - w) A* (uj - w) , 

(18) 

where Z,Z' = {I,Q,U} and a,P = {T, E, B} with 
lo it = 1, uj ue = sin2</> w , oj ub = cos20 w , wq_b = 
cos20 w , ijJqb = — sin20 w and otherwise zero. 

Here we should note that the window functions 
W^ zla o(\w\) are independent of the band-power spec- 
tra and therefore we only need to pre-calculate the in- 
tegrals of the window functions over w in Eq. [17] once 
for evaluating the covariance matrices. Furthermore, 
if the primary beam pattern is Gaussian, the window 
fun ctions can be integr a ted o ut analytically. Follow- 
ing [Hobion~^Maism£ei] (|2002[ ). in Table I2TT1 we provide 
the formulas of integrals for computing the window func- 
tions in Eq. 1181 This is all the formalism required for 
constructing the covariance matrices. 



Empirically, direct evaluation of the full log-likelihood 
function over high-dimensional parameter spaces is un- 
achievable. However, thanks to sophisticated and effi- 
cient numerical algorithms, it becomes possible to find 
the parameter values that maximize the log-likelihood 
function in r elatively few steps, of ord er of Njf. As 
mentioned bv lHobson fc Maisingerl (|2002[ ) and references 
therein, the most efficient numerical algorithm for max- 
imizing the likelihood function is the combination of the 
sparse matrix conjugate- gradient algorithm and Powell's 
directional- set method. Due to the high sparsity of the 
covariance matrix for interferometer data, sparse matrix 
algorithms can dramatically reduce the computational 
time, by a factor of / s 15 compared to the standard dense 
matrix algorithm, where f s is defined by the sparsity frac- 
tion of the covariance matrix. Using Powell's direction- 
set method instead of the standard Newton-Raphson 
method, which requires intensive computations on the 
gradient or curvature of the log-likelihood function, the 
maximising process requires only ~ 3iV 2 function calls by 
the line-minimisation method. For a QUBIC-like obser- 
vation of a single field, with 6 spectral bins in each CMB 
power spectrum, the maximum-likelihood evaluations for 
a total about 4000 visibilities in the I,Q,U maps require 
about 20 hours of CPU time. 

Assuming that the likelihood function near its peak 
a can be well approximated by a Gaussian, the pa- 
rameter confidence intervals can be estimated by tak- 
ing the inverse of the Hessian matrix H(a), which is the 
matrix of second derivatives of the log-likelihood func- 
tion with respect to the parameters, i.e. d 2 H/daidaj. 
The inverse of the Hessian matrix can be regarded as 
the asymptotic covariance matrix of the parameter esti- 
mates. The square roots of the diagonal elements of the 
asymptotic covariance matrix are assumed to be asymp- 
totic standard errors of the parameter estimates, namely 
(Saf) = (H -1 ).;.;. Practically, we perform second dif- 
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Fig. 1. — Simulated interferometric observations. The images shown in panels (a), (b) and (c) are a 30 X 30-degree realization of the 
two-dimensional CMB Stokes fields I, Q and U based on the standard CMB power spectra, with a 64 X 64 pixel grid. All Fourier modes 
higher than the Nyquist frequency are filtered out to avoid aliasing. The images shown in the remaining panels arc simulated Stokes 
visibilities (shown as magnitudes) by a QUBlC-likc observation, assuming a Gaussian primary beam A(x) with beam width cr = 5°, a 
simulated 12-h tin-coverage of single field with 400 close-packed antennas and Gaussian random noise of 0.015/iK per visibility. The map 
units are equivalent thermodynamic temperature in /iK. 

ferences numerically along each parameter direction to 
directly obtain the Hessian matrix and then calculate 
the statistical error estimates on the band-power spec- 
tra. We find this procedure requires only about 30 mins 
of CPU time for ~ 4000 visibilities. 

Fig. [2] shows the resulting maximum-likelihood CMB 
power spectra based on the simulated observations in 
the absence of systematic errors. The recovered power 
spectra are basically consistent with the true underlying 
CMB power spectra within 2-o~. 



tify how significant this systematic error is when com- 
pared with the systematics -free 1-a statistica l error , fol- 
lowing lO'Dea et all (|2007f ) and iMiller et al.l (|2008t ). we 
introduce the tolerance parameters defined by 



AC* 



(20) 



4. ANALYSIS OF POINTING ERRORS 

Using the simulated Stokes visibilities and applying the 
ML analysis described in the previous section, we ob- 
tained estimates of the systematic pointing errors in the 
CMB power spectra. Other systematics, such as beam 
shape errors, gain errors and cross polarization, will be 
presented in a detailed analysis on various instrumental 
systematics in a forthcoming paper. 

The quadrature difference between the recovered power 
spectrum with and without pointing errors is used to es- 
timate the effect of systematic pointing errors. Also the 
pointing errors can potentially change the statistical er- 
ror (which depends on the curvature of the likelihood 
function) for a given experiment. The bias in the i-th 
band-power spectrum Ci and the change in the corre- 
sponding statistical error ~ol are given by 

\2\l/2 



AC; 

Aov 



Ci. 

■^) 2 > 1/2 



(19) 



where Cf rror refers to values obtained in the presence 
of pointing errors and C; is the recovered power spec- 
trum for that same patch of the sky in the absence of 
systematic errors (refers to "ML" in Fig. [2]). To quan- 



where AC; and Ao7 are quadrature differences in Eq. 1191 
We set up a tolerance limit, say 10%, which requires that 
neither a nor j3 exceed 0.1. In our simulations, we con- 
sider two types of pointing errors, which we call uncorre- 
lated and fully- correlated. In the uncorrelated case, the 
pointing errors are assumed to be Gaussian-distributed 
for each antenna (and tend to average out), and in the 
fully-correlated case, the pointing offsets Ax relative to 
the desired observing direction are identical for all the 
antennas and remain fixed on the sky as the sky ro- 
tates. Note this situation is different from the case of 
a misalignment of the entire interferometer by Ax from 
the expected direction, where the observed patch would 
change as the sky ro tates. 

As introduced by iBunnl ([20071 ). we use an error pa- 
rameter p in units of the beam width a to characterize 
the pointing error level, i.e., p = 5/a and p = |Ax|/er re- 
spectively for the uncorrelated and fully-correlated cases, 
where 5 is the root-mean-square (rms) value of Gaussian- 
distributed pointing errors and |Ax| is the amplitude of 
identical pointing offsets. For clear comparison, Fig. [3] 
shows the values of a and j3 in the BB, TB and EB 
power spectra with the identical p = 0.1 for both the 
cases. As illustrated in Fig. [3] for the uncorrelated case, 
the parameter a in each band-power bin for the BB 
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power spectrum is smaller than 0.06 and smaller than 
0.04 and 0.01 for the TB and EB power spectra, respec- 
tively. Also, the parameter a for the BB power spectrum 
increases approximately monotonically with increasing I. 
Moreover, it is apparent that for the BB power spec- 
trum, the values of the parameter (3 are always much 
smaller than those of a over the whole multipolc range, 
except for the lowest band-power bin where (3 roughly 
approaches our criterion of 10% threshold. But for the 
TB and EB power spectra, the values of f3 demonstrate 
that the changes in statistical errors in all band-power 
bins are basically comparable and they have less sensi- 
tivity than a. Using a and /3, and taking the tolerance 
limit of 10%, we put an upper limit on the allowed range 
of the pointing error parameter p, requiring p < 0.1. 

Here a and /3 are estimated by the mean of 50 simu- 
lations, which in principle should be an accurate estima- 
tion of the true systematic contamination. To check the 
convergence of our results, we repeated another indepen- 
dent 50 realizations and found the variations of both a 
and j3 to be less than 5%. In order to clearly illustrate 
the systematic effects, for each simulation we randomly 
generate Gaussian pointing offsets for all the antennas, 
but the Stokes parameters, /(x), Q(x) and ?7(x), and the 
noise in each visibility are fixed during the entire simu- 
lation. In fact, our results are insensitive to different 
realizations of the sky maps and noise and only depend 
on the amplitude of the pointing error. 

To verify our r esults we com pare them with the ana- 
lytical results by iBunrJ (|2007h . In our simulations, we 
assume two visibilities Q and U are measured with the 
sa me an t enna (as is the case of a circular experiment 
in iBunnl (|2007h ). and therefore using the corresponding 
relation, the bias on the i?-mode power spectrum is 

(AC,) 2 = ^ (s(^) 2 (Cf E f + ^Cf E Cf B ) , (21) 

where s 2 is the average of sin 2 2cf> over the antenna pat- 
terns and the factor N is the number of baselines con- 
tributing to each band-power bin since the systematic ef- 
fects in random_pointing offsets will average down as y/N. 
The value of s 2 approximately follows s 2 = 262.7/i 2 
in a QUBIC-like experiment and our simulation shows 
that the number of the baselines in the first two bins is 
about 2000 and about 5000 in the remaining bins. Due 
to pointing errors mostly arising from mixing EE into 
BB, we only consider the leading contribution from Cf E 
and the secondary one from Cf E and Cf B . Using the 
above param eters, the corresponding results are shown 
in Fig. 3(a)| from which we find that the simulation- 
based results are consistent within a factor of 4 with 
the analytical-based results. The simulated observations 
with likelihood-function analysis therefore provide reli- 
able estimates on pointing errors. Furthermore, numer- 
ical sim ulation s illus trate that the analytical-based re- 
sults by [Bunn ( 2003) actually underestimate the biases 
at I > 150, since the analytical calculation can only 
give a first order approximation and thus provide lower 
bounds on pointing errors, which can yield poor esti- 
mates in cases of large pointing errors. Therefore the 
analytical-based results can be considered as an approxi- 
mate estimation of systematic effects. The resulting bias 
that the simulations predict could be larger than those 
found by the analytical calculations. 



In the fully-correlated case, the pointing directions of 
all antennas are offset by the same amount, which means 
they all look at a slightly different patch of the sky from 
the one used to compute Cj in Eq. [501 Recall C, is the 
recovered band-power spectrum in the absence of sys- 
tematic errors. In this case, the systematic errors would 
be equivalent to the cosmic variance and a ~ 1 if in- 
strumental noises were zero. It is worth noticing that 
the fully-correlated case would have no effect on the de- 
termination of cosmological parameters since observing 
a different patch of the statistically isotropic CMB sky 
would not affect the correct estimates of the underlying 
CMB power spectra. The simulations confirm, as ex- 
pected, that the contamination levels in a and /3 in each 
band-power bin are much larger than in the case of uncor- 
rected errors by a factor of ~ 10. In the fully-correlated 
case pointing errors of 10% beam width can affect the 
BB and TB power spectrum measurements at roughly 
the 50% level in the lowest ^-bin. Similar to the uncor- 
rected case, the parameter a is more sensitive to the 
pointing errors than j3. The fully-correlated case how- 
ever is a worst-case scenario and unrealistic. The actual 
pointing offsets should be very close to Gaussian distri- 
butions, which implies that the realistic contamination 
induced by pointing errors would resemble the uncorre- 
cted pointing error forecasts. We conclude that for a 
QUBIC-like observation, the tolerance parameters a and 
P remain below our tolerance limit of 10% for pointing 
uncertainties as large as p = 0.1 for all band powers. 
For the EB and TB power spectra, pointing errors are 
so small that we can entirely neglect such systematic ef- 
fects in all band powers. 

5. CONCLUSIONS 

In this study, we developed a complete simulation 
pipeline to assess systematic errors in measurements of 
the CMB by interferometers. Although we only focus on 
pointing errors at present, any other systematic errors, 
such as beam shape errors, gain errors and cross polar- 
ization, can be evaluated in the same way and we plan to 
study them in a forthcoming paper. The main purpose of 
the present paper is to introduce the maximum likelihood 
estimator for interferometric CMB temperature and po- 
larization data, and to study its sensitivity to incorrect 
modeling assumptions and systematic errors. We choose 
tolerance levels of 10% on a and (3 in this study, which 
are somewhat arbitrary and may be changed at will. For 
interferometers with large numbers of redundant base- 
lines, each independently measuring the Fourier modes 
of the sky, the effects of random systematic errors will 
be highly suppressed since errors will average down in 
visibilities measured over many baselines. 

For a QUBIC-like interferometer, we find that in most 
cases the most stringent constraints on the allowable 
range of pointing errors are obtained from the require- 
ment on the bias in power spectrum rather than from 
the changes in statistical errors. When the Gaussian- 
distributed pointing errors are controlled with a preci- 
sion of S ss 0.7°, our simulation shows that the measured 
-B-modes in the multipole range 28 < I < 384 can not be 
contaminated at the 10% level in the statistical uncer- 
tainty (ae) units, but the change in the statistical error 
at 28 < £ < 88 could exceed 10% of the statistical un- 
certainty. Our results are consistent with the analytical 
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Fig. 2.— The CMB power spectra TT, EE, BB, TE, TB, EB (red) recovered by maximizing the likelihood function from a mock QUBIC- 
like observation in the absence of systematic errors (see details in text). The fiat band powers are estimated in 6 bins with bin- widths 
Au = 9.5 and with 1-a (red hatched) statistical uncertainties. The input CMB power spectra are shown in black and the specific realizations 
of such input power spectra are shown in green. 





200 300 
Multipole 1 



(a) BB power spectrum 



o.ooi 




o.ooi 




Multipole 1 

(b) TB power spectrum 



Multipole 1 

(c) EB power spectrum 



400 



Fig. 3. — Assessment of systematic pointing error for BB, TB and EB power spectra. In the units of systematics-free 1-cr statistical 
uncertainties, the solid lines are for the induced bias in the power spectra and dotted for the changes in statistical uncertainties. Green 
triangles denote the results in the case of the independent Gaussian-distributed pointing errors with the di s persio n 8 = 0.1a for each 
antenna, which are estimated by averaging over 50 simulations. The corresponding analytical result by [Burin (200719 ' s shown in dashed- 
black line. Blue points correspond to all the antennas having the identical pointing offsets with the amplitude |Ax| = O.ler and the results 
are obtained by averaging over 20 random directions relative to the target direction. 
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estimations by iBunnl ()2007l ). within a factor of 4. 

Moreover, the TB and EB power spectra are unique 
"smoking gun" signals for new physics in principle. For 
imaging experiments, pointing errors have to be con- 
trolled to the sub-arcminute level to avoid spurious sig- 
nals. However, our simulation clearly shows that the 
impact of the pointing errors on TB and EB estimates 
are negligibly small compared with their statistical un- 
certainties. Therefore systematic pointing errors in in- 
terferometers will not severely degrade B-mode science. 
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